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1. Introduction 

The development of classical phase space methods to describe quantum systems whose 
kinematics is governed by Cartesian variables has its origin in two important and 
independent ideas. The first, due to Weyl [l] is the setting up of an unambiguous 
rule or convention that maps each (real) classical dynamical variable, a classical phase- 
space function, into a corresponding (hermitian) operator for the quantum system in a 
linear manner. The second, due to Wigner [2] , is the definition of a (real) phase-space 
distribution function representing each (pure or mixed) quantum state in a complete 
manner. Later work [Hj clarified that these two rules or definitions are exactly inverses of 
one another, as a result of which quantum mechanical expectation values can be written 
in a classical-looking form as phase-space integrals. 

In an interesting paper developing the analogies between classical and quantum 
mechanics, Dirac [3] discussed the general problem of expressing a quantum dynamical 
variable, an operator, as a function of the basic complete and irreducible set of operators 
of the quantum system. The latter forms a non-commutative set, leading to the concept 
of ordering rules in forming functions of them. In this context, Dirac used a description 
of an operator by a collection of mixed matrix elements, the rows and columns labelled 
by two different orthonormal bases for the Hilbert space of the system. 

Phase space distributions have played a significant role in optics, particularly in 
unifying radiometry and radiative transfer with the theory of partial coherence. These 
activities were triggered by the pioneering work of Walther on radiometry and that of 
Wolf on radiative transfer. Walther introduced two definitions for the radiance function. 
The first definition [3] is analogous to the Wigner distribution and, indeed, the reader's 
attention was drawn by the author to this fact. Walther's second definition which 
has since been used in hundreds of radiometry papers may be seen, in retrospect, to 
be analogous to the Dirac inspired view of phase space distributions developed here. 
Interestingly, this remark applies equally well to Wolf's expression for the specific 
intensity , having been inspired by the second definition of Walther's. 

The purpose of the present paper is to show that one can arrive at the Weyl- 
Wigner formalism and results starting from Dirac's ideas and following a physically 
well-motivated and at the same time an extremely elementary and transparent series of 
steps. This illuminates the use of phase-space language for quantum mechanics from 
an alternative perspective. The properties of Wigner distributions have been studied 
in considerable detail by many authors. It has been shown for instance, that if one 
asks for a phase-space description of a quantum state obeying a small number of very 
reasonable conditions, then the Wigner distribution is the unique answer. In important 
work by several authors [H] the Weyl- Wigner formalism has been shown to be one of 
several different possibilities all based on the use of phase-space methods. In the light of 
this, an alternative line of argument leading directly and very easily to the same answer 
may be of interest. 

Among the important properties of the Wigner distribution corresponding to a 
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given quantum state is that one recovers the two quantum probability distributions 
in position or in momentum space by integrating over the momentum or the position 
variable respectively, in the Wigner distribution. This is called the recovery of the 
marginals. In our treatment this too becomes immediately evident upon inspection 
and completely transparent, no calculations being needed at all. The value of 
expressing things in this way is that it may suggest possible generalizations to other 
quantum-mechanical situations where the basic variables and commutation relations 
are not of the Cartesian Heisenberg type. We have in mind for instance the case of 
finite-dimensional quantum systems jHj, quantum mechanics on Lie groups jTU], non- 
commutative geometric quantum schemes, deformed or g-quantum kinematics etc. In 
all these cases the presentation of the very well known and familiar case in a very concise 
and transparent manner may help suggest interesting modifications to be studied. 

A brief summary of the present work is as follows. After a quick review of the 
familiar phase space descriptions of operators in quantum mechanics in Section 2, in 
Section 3 we examine the question as to how the trace of products of two operators can 
be expressed as phase-space integrals in terms of their phase space representatives in 
such a way that the inherent symmetry of the trace operation is manifest. We show that 
this can be done at the expense of introducing a kernel. We investigate the properties 
of this kernel in detail and show how by defining a square root of the kernel in a special 
way one is naturally lead to the definition of Wigner distribution and associated phase 
point operators. Some consequences of the Wigner- Weyl correspondence are examined 
in Section 4, followed by a discussion on the recovery of the marginal distributions in 
Section 5. Section 6 contains our concluding remarks and further outlook. 

2. Elementary Phase-Space Descriptions of Operators 

We consider a one-dimensional quantum system whose basic operators are a hermitian 
Cartesian pair q, p obeying the Heisenberg commutation relation: 



For corresponding classical phase-space variables, as well as particular classical values 
or quantum eigenvalues, we use q, q', q", ■ ■ ■ ,p,p',p", ■ ■ ■. The continuum-normalized 
eigenstates \q), \p) of q and p obey as usual: 



[q, p\ = ih. 



(1) 



(q\q') =5(q-q'), (p\p) = 8 (p - p') , (q\p) = (2nh) 1/2 exp (iqp/ K) . (2) 



From the completeness relations: 




(3) 



— oo 



— oo 



we obtain immediately the two operator statements: 



£ (q- q) = \q)(q\, s (p-p) = \p)(p\. 



(4) 
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Consider now a quantum-mechanical operator A. It can certainly be completely 
described by its position-space matrix elements {q'\A\q) which constitute in general a 
non-local kernel (In case A is unitary, (q'\A\q) is the overlap between eigenstates of the 
"old" position operator q and a "new" one: q 1 = Aq'A' 1 , and then the kernel of A 
would be the exponential of (i times) the analogue of the classical generating function 
of a canonical transformation of the type "q — Q" PI]). The kernel corresponding to A^ 
is: 

(q'\$\q) = (q\A\q'y. (5) 

To move towards a description of A at a classical phase-space level it is natural to 
consider, in the spirit of Dirac [H], the mixed matrix element (g|A|;p) which, regarded 
as a function of the phase-space variables q and p, certainly also describes A completely. 
For later convenience we include a non-vanishing plane wave factor and define the "left" 
phase-space representative of A as the function: 

A (q,p) = (q\A\p)(p\q) = Tr I A\p) (p\q) (q\ 



= Tr{A5(p-p)5(q-q)j (6) 
= (2nh) (q\A\p) exp (—iqp/h) . 

Here, of course, (p\q) is the kernel of the unitary operator corresponding to Fourier 
transformation, which interchanges q and p. It is interesting to note that in Dirac's 
treatment ^T] Ai(q,p) is regarded essentially as the ratio (q\A\p) / (q\p) , which 
determines the form of A as a function of q and p in standard ordered form, i.e. 
q to the left of p. 

Even if A is hermitian, Ai (q,p) is in general complex. However we do have, as is 
particularly obvious from the second line in Eq.©: 



dpAi (q,p) = (q\A\q), / dqA x (q,p) = (p\A\p), 
Tr iA> = / / dqdpAi(q,p) . 
As an alternative to the above, the "right" phase-space representative of A is given 

by: 



A (q,p) = (p\A\q}(q\p) = Tr [A\q) (q\p) (p\ 

= Tr\A8(q-q)8(p-p)\ (8) 



= (2nh) 1 (p\A\q) exp (iqp/h) . 

This is related to expressing A in anti-standard form, p to the left of q, and we have 
again: 

dpA r (q,p) = (q\A\q), \ dqA r (q,p) = (p\A\p), 

(9) 



Tr |a| = / / dqdpA r (q,p) 
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Thus we have two equally elementary phase-space descriptions of the operator A on the 
same footing, with the roles of coordinate and momentum interchanged to go from one 
to the other. As noted above, even in the hermitian case in general both A\ (q, p) and 
A r . (q,p) are complex. More generally under hermitian conjugation we have: 

B=:A^B r (q,p) = A l (q,p)*, (10) 

so in the hermitian case we have: 

fi = A=i>A r (q,p)=A l (q,py. (11) 

We can now ask if we can pass in a natural way to a third phase-space description 
of A standing exactly "midway" between Ai (q,p) and A r (q,p), thus treating q and p 
symmetrically. This is achieved in the next Section. 



3. Operator Product Traces and Passage to the Weyl- Wigner Description 

Consider two generally non-commuting operators A and B. The trace of their product 
is symmetric under their interchange and can be expressed in two ways using classical 
phase space§: 

Tr{A§} = J J dqdp(q\A\p)(p\B\q) 

= 2nhJ J dqdpA { (q,p) B r (q,p) (12) 

= 2-Kh J J dqdpA r (q, p) Bi (q, p) . 

The last line follows from the previous one by symmetry under interchange of A and B. 
However in each of these two phase-space integrals the manifest symmetry in A and 
B is lacking. One can ask if such symmetry can be restored while continuing to work 
with phase-space quantities. Towards this end we begin by first expressing (AB)i (q,p) 
entirely in terms of Ai (q',p') and B\ (q",p")' 



(ABj^p) =(q\AB\p)(p\q) = j J dq'dp' {q\A\p')(p'\q'){q'\B\p)(p\q) 

= J J dq'dp'A l (q, p') K x (q, p'; q> , p) B t (q> , p) , (13) 

Kl(qM,p) = (^h) 2 (q\p')(p'\q')(q'\p)(p\q)=e W {i(q-q')(p'-p)/h}. 

the first line in the definition of Ki following from: (q\A\p') = Ai(q,p') /(p'\q) and 
from: l/{p'\q) = 2nh(q\p'). The non-local convolution involved in expressing {AB) l in 
terms of Ai and B\ is an indication already of the general situation since we are dealing 
with mixed matrix elements; it is a forerunner of the Moyal or "star" product when 
the transition to the Weyl- Wigner description of operators is completed. We may also 
note that, aside from the continuum normalization of the q and p eigenvectors, the 

§ Thus, as is well known, from the point of view of trace calculations the standard and anti-standard 
orderings are dual to one another. 
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kernel K\ is a four- vertex Bargmann invariant ^2]- Hence its phase, which is the area 
of the phase-space rectangle with vertices (q,p) , (q,p') , (q',p) and (q',p r ) is a geometric 
phase [IB]- Now combining Eq. ljlHj) with Eq.((ZJ) and relabelling some variables for 
convenience, we get Tr{AB} entirely in terms of left representatives: 



Tr {Ifi} = I I I I dqdpdq'dp'A l (q,p) K x {q,p; q',p') B x {q',p') . (14) 



The kernel K\ {q,P',q',p') is explicitly symmetric under: (q,p) < ► (q',p'), so we have 

a classical phase-space expression for Tr{AB} manifestly symmetric in A and B, but 
at the cost of a kernel. In addition to symmetry, this kernel possesses two important 
properties: it is invariant under phase-space translations as it depends only on the 
differences q — q', p — p'\ and it satisfies the "marginals" equations: 



J dp'Ki (q,p; q',p') = 2ttM (q - q') , 
J dq'Ki (q, p; q', p') = 2vrM (p - p') . 



(15) 



The most natural question is to ask if this kernel can in some sense be "transformed 
away" while maintaining manifest symmetry in A and B. This can be done if we can 
express it as the "square" or the convolution of some more elementary kernel, say in the 
form: 



Ki (q, p; q\ P ) = J J d( fdp"t (q", p"; q, p) £ (q", p"; q\ p') . (16) 

From the known properties of Ki(q,p;q',p') we can demand that £ (q , p; q' , p') too 
be a symmetric function of its (pairs of) arguments; be invariant under phase-space 
translations and so depend only on the differences q — q', p — p'\ and possess the 
"marginals" property: 

/ dp'£, (q,p; q',p') = V2%h5 (q - <f) , 

J (17) 
/ dq'£ (q, p; q\ p') = V27ph5 (p - p') . 

Easy calculation shows that the expression: 

t (q,p; q',p') = ^ ex P i 2i (? - <?') (p - p') l n ) ( 18 ) 

obeys all the conditions imposed above||. If we use this in Eq. (jl4)l and associate one 
factor of £ each with A\ and B\ we arrive at the simpler expression: 

Tr{AB} = ^_ J J dqdpA(q,p)B(q,p), (19) 

where A (q, p) arises from A\ (q, p) via: 

A (q, p) = J J dq'dp'i (q, p; q', p') M (q', p') (20) 

| It is interesting that while as a function of four phase-space variables £ is the pointwise square of Ki, 
as a kernel it is the square root of Ki. 
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and similarly for B(q,p)%. We have thus achieved, by a two-step procedure, our 
objective of expressing Tr{AB} as a manifestly symmetric classical phase-space integral, 
with one phase-space function each representing A and B and with no additional kernel. 
One can now see that in the case A* = A, B' — B, since Tr{AB} is real and A and B 
can be chosen independently, A (q, p) and B (q, p) must be individually real: 

A = A^A(q,p)=A(q,p)*. (21) 

Expression ()20|) for A (q, p) is indeed the Weyl- Wigner representative of A in 
phase-space form. With elementary manipulations we can express it in the equivalent 
forms [Tij : 

A(q,p) =2 J J dq'dp'A l {q 1 , p') exp {2i {q - q') (p - p') /h} (22) 
dq'(q - ^q'\A\q + -q') exp {ipq'/h} (23) 



or: 



where: 



A{q,p) = 2nhTr {AW(q,p)} , (24) 

W(q,p) =7^1 J dq'\q+^q')(q-^q'\exp{ipq'/h} (25) 

dq'\q + -q')(q + ^q'\p)(p\q ~ ^q')(q - ^q'\- 

Representing this "phase-point operator" as: 

VV(q,p) = J dq'dq"\q')(q'\W(q,p) \q")(q"\, (26) 

the matrix elements are given by: 

(q'\W (q,p) \q") = \d (q' + q" - 2q) exp {ip (q> - q") /%} . (27) 
Tin 

It is evident from the first line in Eq. (j25)l that W {q,p) is hermitian as well as that: 

Tr{vV(g,p)} = ^. (28) 

The Weyl correspondence makes use (see below) of the "Weyl operators" 
exp{i(xp — kq) /h} that are labelled by the phase-space points x and k. Then, one 
can prove that: 

W(q,p) = j J (2^)2 exp ^ _ ^ xp ~ kq)}exp{i(xp- kq) /h} , (29) 

i.e. that the phase-point operators are the symplectic Fourier transforms of the Weyl 
operators. Indeed, a straightforward calculation shows that: 

(q'\ exp {i (xp — kq) /h} \q") = 5 (x + q — q") exp < — ik ^ — - i (30) 



2h 

and using this result to evaluate the matrix elements of both sides of Eq. fl29|) one obtains 
back Eq.(j2ZD. 

We look at some more properties of the operators W (q, p) in the next Section. 
% The factor in Ea. (|19f) is chosen so that A and A (q,p) have the same physical dimension. 
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4. Some Consequences of the Weyl- Wigner Correspondence 

It is illuminating to compute the phase-space functions representing the operators 
<5 (<? — <?) 8 (p — P) i W (<LP) an d <5 ip — V) & (<? ~~ <?) occurring in Eqs.tjHJ), © and (f2l)j) 
under the Weyl- Wigner correspondence. This helps us also understand the structure 
of the kernel £ (q,P',q'p') better. It is well known that the Weyl correspondence takes 
elementary classical exponentials to elementary operator exponentials according to: 

exp {i (aq — rp)} — > exp {i (aq — rp)} (31) 

and then extends this to general functions A(q,p) by linearity. For the operators 
mentioned above we find using Eqs.flBJ) and (J2Uj) : 

A = 8(q-q')8(p-p') ^Ai{q,p) = (q - q') 8 (p - p') 

Znn 

^A{q,p) = \ exp {2i (q - q') (p - p') /%} , (32) 

A = 5(p- p') 8{q-q')^A (q,p) = \ exp {-2i (q - q') (p - p') /H} , (33) 

Tin 

the second result following from the first by hermitian conjugation, and we see that, up 
to a factor, the latter is just the kernel £ (q,p; q',p') of Eq. (fTH|) . For W {q',p') we find 
from Eq. (j25} : 

A = W (q', p') Ai (q, p) = — exp {-2i (q - q') (p - p') /%} 

^A(q,p)=5(q-q')5(p-p > ). (34) 

At first sight, it is not easy to recognize that the operators W(q',p') stand "midway" 
between 5 (q — q') 5 (p — p') and 5 (p — p') 5 (q — q'), in the sense of treating q 
and p symmetrically, or that the classical representative 5 (q — q') 5 (p — p') of 
W (</,£>') also stands "midway" as a phase-space function between the two functions 
(l/Tih) exp {±2i (q — q') (p — p') /h}. However this is actually so, as can be appreciated 
by looking at the Fourier transforms with respect to q' and p' . Classically we have: 

— p J J dadrexp {i[a (q — q') — r {p — p')]} x {exp (±ihcrr/2) or 1} 

= — exp {±2i (q — q') {p — p') /h} or 5 {q — q') 5 (p — p') (35) 

TXfi 

and the Weyl map then preserves these "relative positions" among the corresponding 
operators. 

Two other known properties of the W(q,p) are immediately read off from the 
relations assembled above, with no need for any calculations. Since 5 (q — q') 5 (p — p') 
is real, W(q',p') is hermitian, as we have already noted, and from Eqs.(|THJ) and (|3~4"j) 
they are trace orthonormal in the continuum sense: 

Tr | W (q, p) W (q, p) j = ^ J J dq"dp"5 (q - q") 5 (p - p") 5 {q - q") 5 (p - p") 

= 2^h 6{q ~ q>) 6 {p ~ pl) • (36) 



(2 



Phase-space descriptions of operators and the Wigner distribution I 



9 



Thus the inverse of Eg . (J24j) reads: 

A = [ / dqdpA(q,p)W(q,p). (37) 



It is known [15] that the operators W (q, p) obey the following interesting anti- 
commutation relations with q and p: 

' ( Q, VV (q, p) } = qW (q, p) , - \p, VV (q, p) } = pW (q, p) . (38) 



2 I 11 y171 ' J * v "' x y ' 2 
These are operator versions of corresponding immediately obvious classical relations 
stemming from: q5 (q — q') = q'5 (q — q'). If we set q = p = in Eq. (j38j) we see that 
W (0, 0) anticommutes with both q and p: 

g VV (0, 0) = -VV (0, 0) q, p W (0, 0) = -VV (0, 0) p. (39) 

In turn this means that W (0, 0) 2 commutes with both q and p, so it must be a multiple 
of the unit operator. We can easily convince ourselves that W (0, 0) is nonzero, so 
we must be dealing here with a nonzero multiple, which means that W (0, 0) has an 
inverse. The existence of the inverse together with the first of Eqs. (j39j) leads at once 
to: VV (0, 0) q VV (0, 0) _1 = — q and similarly for p, so the upshot is that VV (0, 0) is a 
multiple of the parity operator. Now, e.g. in the {\q}} basis, the matrix elements of the 
parity operator P are given by (cfr.Eq.(j2J)): 

(q\P\q') = (q\-q')=5(q + q') (40) 

and then, in this basis: 

Tr{P} = J dq(q\P\q) = J dqS (2q) = \. (41) 



Taking then traces and using Eq. (j28|) fixes the proportionality factor and we find: 

VV(0, 0) = —P. (42) 

Summarizing, one finds the Weyl correspondence: 

A (q, p)=6(q)6(p)^A = W (0, 0) = \p. (43) 

Remark. We would like to stress that what we mean here by "trace" of an operator 
is defined as the sum (or integral) of the diagonal matrix elements of the operator with 
respect to a given basis, and that the only bases that we are using here are the {\q}} 
and {\p)} bases, in each of which, e.g. and consistently: Tr{P} = 1/2. This does not 
mean that we are claiming that P is a trace-class operator, which would imply much 
stronger requirements that the parity operator is not likely to meet. 

One can carry the previous analysis one step further. The unitary operator: 

U (q,p) ='■ exp {i \pq — qp\ /K} (44) 
acts as a phase-space displacement operator. Indeed: 

U(q,p)qU{q,p)* = q-q, U (q,p)pU (q,p) j = p - p (45) 
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and hence, for any (analytic) operator O = O (q, p): 

U (q,p)0(q,p)U (g,p) f = 0(q-q,p-p). (46) 
Taking then: O (q,p) = W (0, 0) from Eq. (j29j) one obtains immediately: 

U (q, p) VV (0, 0) U (q, p? = VV (q, p) (47) 
and hence, from Eq. (|43j) : 

W(q,p) = \u(q,p)PU(q,p^ . (48) 

Thus the phase-point operators are just the parity operation with respect to general 
phase-space points. This leads to the at first sight unexpected operator property: 

W{q,pf = -^l, (49) 

which means that the eigenvalues of VV (q,p) are ±l/nh. This will be used in the next 
Section. 

The purpose of this discussion was to show that this otherwise rather unexpected 
fact is an immediate consequence of the relevant operator relations given above. 

5. Recovery of Marginal Distributions 

We have mentioned in Section 1 that an important property of the Wigner distribution is 
that upon partial integration over either p or q the complementary quantum-mechanical 
probability distribution emerges. In this Section we show how this happens practically 
automatically or manifestly if we use the relations given in Section 3. 

The general Weyl association is as given in Eqs. (j2U|) . (|25|). For the density operator 
p representing a (pure or mixed) quantum state, we use a different normalization and 
define the Wigner distribution by: 

p(q,p)=Tr{pW(q,p)}, (50) 
so that the general expression for the expectation value of A in state p has the form: 

Tr = J J dqdpA(q,p) p(q,p) . (51) 

We notice at this point that Eq.(|49|) concerning the spectrum of W(q,p) has the 
following implication: 

\p(q,p)\<\- (52) 

This known property of Wigner distributions is usually obtained by using the 
Cauchy-Schwartz inequality, so it is interesting to see it emerging here by a much more 
elementary argument. Now we turn to the marginals. 

For a general operator A we have seen in Eqs. (jTj) and (JOJ) that: 

dp\ r (q,p) = (q\A\q), / dqA hr (q,p) = (p\A\p). (53) 
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At the next step, for the kernel £ (q,p; q',p') of Eq.ffTp]) we have the properties in Eq. (fT7|) . 
Combining the above two pairs of equations with the passage (J2U|) from Ai (q, p) to the 
Weyl representative A (q,p) of A, it is immediately seen that: 

dpA(q,p) = 2nh(q\A\q), J dqA(q,p) = 2nh(p\A\p). (54) 

Essentially no work has to be done to get these results. In the case of the density 
operator p we omit the factor 27th, so we recover the marginal probability distributions 

as: 

J dpp(q,p) = (q\p\q), J dqp(q,p) = (p\p\p). (55) 

We can also obtain operator forms of such statements for W (q,p). Knowing that 
the Weyl map is linear, that W {q',p') corresponds to 5 (q — q') 5 (p — p'), and that the 
operators corresponding to 5 (q — q') and 5 (p — p') are S(q — q f ) and S(p — p'), we 
have: 

J ' dp'5{q-q')5{p-p') = 5{q-q') J dp'W (q',p') = 5 (q- q') = \q')(q'\, (56) 

J dq'5 (q - q') 8 {p - p') = 5 {p - p') & J dq'W(q',p') = 5(p-p') = \p')(p'\. (57) 

The purpose of showing the recovery of marginals in this manner is again to emphasize 
how elementary and transparent the derivations are. 



6. Concluding Remarks 

To conclude, we have shown how by expressing the trace of product of two operators 
in terms of their phase space representatives and some fairly elementary steps one is 
naturally led to the concept of the Wigner distribution. Crucial to this construction 
is a kernel with the structure of a Bargmann invariant and its square root endowed 
with certain desirable properties. A noteworthy feature of the approach developed 
here is its economy - no auxiliary constructs are required at all and above all the 
facility with which it lends itself to application to non Cartesian situations such as 
finite state quantum systems which are of particular relevance to quantum computation 
and quantum information processing. That this is indeed the case will be demonstrated 
in a companion paper. 
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